A Q-Ising model application for linear-time image segmentation 
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A computational method is presented which efficiently segments digital grayscale images by di- 
rectly applying the Q-state Ising (or Potts) model. Since the Potts model was first proposed in 1952, 
physicists have studied lattice models to gain deep insights into magnetism and other disordered 
systems. For some time, researchers have realized that digital images may be modeled in much 
the same way as these physical systems (i.e., as a square lattice of numerical values). A major 
drawback in using Potts model methods for image segmentation is that, with conventional methods, 
it processes in exponential time. Advances have been made via certain approximations to reduce 
the segmentation process to power-law time. However, in many applications (such as for sonar 
imagery), real-time processing requires much greater efficiency. This article contains a description 
of an energy minimization technique that applies four Potts (Q-Ising) models directly to the image 
and processes in linear time. The result is analogous to partitioning the system into regions of four 
classes of magnetism. This direct Potts segmentation technique is demonstrated on photographic, 
medical, and acoustic images. 



I. INTRODUCTION 

Segmenting digital images into regions of distinct types 
has appHcations in a great many fields, e.g. medical 
imaging [1-4], surveillance by synthetic aperture radar 
[il and satellite, and underwater acoustic imaging (.e.g. 
see Ref.lli), to name a few. In general, image segmenta- 
tion is performed by classifying image regions by color, 
intensity, and texture. Of course, only the latter two are 
considered in grayscale segmentation, which is the sub- 
ject we consider in this article. While classification by 
intensity is a straightforward assessment of the bright- 
ness/darkness of an image pixel or group of pixels (as 
with histogram segmentation methods [l|), texture clas- 
sification is much more complex [l], [8[ . Although it is dif- 
ficult to define image texture precisely, we can say that it 
is the spatial relationship of the intensities (i.e. "grain- 
iness") of an image region. Clearly, there are indeed a 
large number of possible texture types in digital images. 

A grayscale digital image may be represented as a ma- 
trix of numerical values, as in Fig. [TJ which indicate the 
intensity or brightness (gray level) of the corresponding 
image pixel. Techniques for image segmentation (such as 
thresholding and histogram methods) that focus solely 
on intensity are the most computationally efficient since 
they generally require just one or two passes through the 
intensity matrix. Identifying regions of different texture, 
however, tend to be relatively expensive in terms of com- 
putation time Q. Here I describe four types of mag- 
netism using the Ising model [oj and show how its gener- 
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alization, the Q-Ising (or Potts) model [10], can be used 
analogously to segment grayscale images into four cate- 
gories [11]. The common approach (e.g. see (sl, [6l,fT2l,[T3|) 
to segmentation using the Potts model applies the Potts 
model to the categorizing labels, i.e., the segmented im- 
age. (A similar approach is also used for more general 
data clustering [3, [l5|.) The optimal segmentation is 
often found using a Monte Carlo approach, which pro- 
cesses in exponential time (time for computer processing 
increases exponentially with the number of pixels con- 
tained in the digital image [12]). Use of the Bethe ap- 
proximation [12, 16] allows an optimal segmentation to 
be obtained in power-law time. However, the direct Potts 
segmentation method presented in this article, which 
utilizes both intensity and texture, processes in linear 
time-among the fastest methods currently used. (Other 
unrelated linear-time or nearl y li near-time segmentation 
methods are described in |17H19|. though, these contain 
undesirable feature for many applications. For example, 
graph partitioning methods [18, 19] cannot relate disjoint 
sets with equivalent textures and require training sets.) 
This computational speedup arises from applying the 
Potts model directly to the original image itself rather 
than to the segmented representation. The speedup is 
important for high-data- volume/high-resolution applica- 
tions where time is constrained. 



II. ISING MODEL 

In 1925, Ernst Ising proposed the model H, [20|, [III, 
which now bears his name, as a means to study phase 
transitions in magnetism. Ising solved the original Ising 
model (only two states-namely 1, or spin up, and -1, 
or spin down) in one dimension. However, the one- 
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FIG. 1: On left, a grayscale image varies from black to white 
in going from top to bottom. The image is represented as a 
matrix of intensity values as illustrated on the right. 



dimensional Ising model did not exhibit a phase transi- 
tion, and the two-dimensional Isin g m odel was not solved 
until 1944. At that time Onsager [22| demonstrated that 
the two-dimensional Ising model did, in fact, exhibit an 
order /disorder phase transition. The Ising model and re- 
lated models have been widely used to study magnetism 
and other order/disorder phenomena in physical systems. 
The two-dimensional model is described here, and it will 
be used to define four classes of magnetism. The image 
texture analogs will be presented in a subsequent section. 

In the Ising model, given a two-dimensional lattice of 
magnetic spins oriented either up or down as in Fig. [2] 
(numerically represented as a matrix of I's and — I's), 
the total energy E associated with the lattice is 



FIG. 2: A spin matrix is shown on the left with arrows in- 
dicating elements whose magnetic states are either spin-up 
and spin-down. The matrix on the right represents the same 
matrix with the arrows exchanged for I's and (-l)'s. In the 
absence of an external magnetic field (H = 0), a ferromag- 
net's spins will spontaneously align as depicted in the top half 
of the matrices, while the spins for antiferromagnets will op- 
positely align (alternate in opposite directions) as depicted in 
the bottom half of the matrices. 
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E = -jiH^Si - J ^ SiSj, (1) 

i ij, nn 

where i represents a magnetic element, Si is the value 
of the magnetic spin of i (1 or —1), /i is the magnetic 
moment of the element, H is the external magnetic field, 
and J is the interaction coupling constant. The sum over 
all nearest-neighbors pairs ij is represented by 
Nearest neighbors are considered to be the elements di- 
rectly above, below, left, and right of a given element. 
If no coupling interaction is present (J = 0) then ja > 
creates an energetically favorable condition for the spins 
to align in the direction of the external field whereas 
/i < leads to spins aligned in the direction opposite of 
H. (See Fig. [3) 

III. FOUR CLASSES OF MAGNETISM 

In common language, magnetism usually refers to the 
behavior of magnets that attract iron (known as fer- 
romagnets). Other materials, however, reveal different 
magnetic behavior. These behaviors may be classified 
into the following four groups: 

1. Ferromagnetism is exhibited by materials for 
which, even in the absence of an external magnetic 



FIG. 3: Magnetic spins are aligned with (left, top) and against 
(left, bottom) the magnetic field H. These regions are asso- 
ciated with paramagnetism and diamagnetism, respectively. 

field, a net magnetization is present with north and 
south poles. A ferromagnet will both attract and 
repel another ferromagnet depending on the align- 
ment of their poles. 

2. Antiferromagnetism is displayed by materials 
which do not form magnetic poles even with a mod- 
est external magnetic field present. An antiferro- 
magnet will neither attract nor repel another anti- 
ferromagnet. Likewise an antiferromagnet will have 
no magnetic interaction with a ferromagnet. 

3. Paramagnetism shows no spin ordering in the ab- 
sence of an external field, although the paramag- 
net's spins align in the same direction of an exter- 
nal field when present. A paramagnet will neither 
attract nor repel another paramagnet, but will be 
attracted to a ferromagnet 

4. Diamagnetism also shows no ordering in the ab- 
sence of an external field; however, the diamagnet's 
spins will align in the opposite direction of any ex- 
ternal field. A diamagnet will neither attract nor 
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repel another diamagnet and will be repelled by a 
ferromagnet. 



spins that may or may not belong to Aa. Then the total 
energy E is simply the sum over all regions, 



IV. Q-ISING MODEL 

This section describes the Q-Ising model, also known 
as the Potts model, which will be used to analyze 
grayscale images where the pixel intensities are consid- 
ered in place of the magnetic "spin" values. The Potts 
model generalizes the (2-state) Ising model into a Q-state 
model by replacing the I's and -I's with the numbers 
0, 1,2,3, ...,(5 — 1 (the pixel intensities in a grayscale 
image). According to this model, the total energy asso- 
ciated with the system is 



(7) 



E = -iiH Y,{si -st)-jJ2 



(2) 



where St is the spin threshold (separating light from dark 
with < St < Q — I) and the kronecker delta function 
^ (1 if the elements are equal, otherwise) replaces the 
summand for the second summation in Eq. [TJ Now lets 
suppose that a magnetic spin matrix represents the in- 
tensities in a grayscale image and that a uniform external 
magnetic field, H = 1, is applied to the system. We wish 
to find a set of pairs (/i^, J^) for Eq. [2] which character- 
izes the regions a in the spin matrix. Then each spin in 
the matrix will be associated with a particular parameter 
pair. The "optimal characterization" is taken to be the 
partitioning of the spin matrix in such a way that the to- 
tal energy of the matrix is minimized. Minimizing Eq. [2] 
requires the set of four parameter pairs (/i > 0, J > 0), 
(// > 0, J < 0), (/i < 0, J > 0), and (/i < 0, J < 0). For 
convenience, and without loss of generality, we simply use 

(/il = 1, Ji = 1), (/i2 = 1,-^2 = -1), (M3 -1,-^3 1), 

and (/i4 = — 1, J4 = —1) for what we will call regions 
a = 1,2,3,4 respectively. The energies for the regions 
can be written as 



Ei = -Y,isi-St)- (3) 

ieAl ieAljenn 

E2 = -Y,{si- St) + •^(«'' (4) 

ieA2 ieA2jenn 
ieA3 ieASjenn 



^4=^(Si-St)+ Y <^(«i'«i)' (6) 

ieA4 ieA4,jenn 

where Aa is the set of all spins i labeled with category 
a. The j's in the interaction terms are nearest-neighbor 



where a = 1,2,3,4 are the region segments assigned to 
categories 1, 2, 3, and 4, respectively. These energies cor- 
respond to regions which are 1) paramagnetic and fer- 
romagnetic, 2) paramagnetic and antiferromagnetic, 3) 
diamagnetic and ferromagnetic, and 4) diamagnetic and 
antiferromagnetic. Minimizing the total energy (Eq. [7]) 
effectively partitions the spin matrix into regions corre- 
sponding to these four classes of magnetism. Reverting 
back to the grayscale image, these categories correspond 
to image regions which are 1) bright and smooth, 2) 
bright and "grainy" , 3) dark and smooth, and 4) dark and 
"grainy" . The direct Potts segmentation process requires 
the matrix to be separated into regions of these four types 
in such a way as to minimize the total energy. An exam- 
ple image with four distinct intensity and texture regions 
is depicted in Fig. and the result from the direct Potts 
segmentation is shown in Fig. |4)3. The parameters used 
in this example are Q = 16 and St = 7.5, and the seg- 
mented image is shown with four shades of gray repre- 
senting the four texture categories. The segmented image 
is filtered using the ImageMagick[28j despeckle algorithm 
and is shown in Fig.Hl:. For comparison, a segmented im- 
age using a simple four-level histogram method is shown 
in Fig. mi. The four-level histogram method performs 
poorly in distinguishing the image texture regions. 



V. ALGORITHM 

To segment a digital grayscale image, we first trans- 
form the image into a matrix of intensity values Sij 
as illustrated in Fig. [H Since it is common for 8- 
bit digital grayscale images to contain 256 gray levels 
(0 < Si J < 255), reducing the gray levels to 16 < Q < 64 
usually helps in identifying visual textures. Next we must 
specify an intensity threshold 5^, which separates bright 
pixels from dark pixels. For each of the examples in this 
article, the threshold is simply St = (Q — 1)/2, or halfway 
between maximum and minimum intensity. (If texture 
identification at a certain scale is desired, one may rescale 
the image.) Then to perform the energy minimization de- 
scribed above, we follow these two steps for each of the 
image pixels (i, j): 

1. If Si J > St then (i, j) is assigned either texture cat- 
egory 1 or category 2. Otherwise, (i, j) is assigned 
texture category 3 or 4. 

2. Check the nearest-neighbor values (i.e., the pixels 
directly above, below, left, and right) of pixel (i^j). 
If one or more nearest-neighbor values are equal to 
Si J, then pixel (z, j) is assigned texture category 
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(c) (d) 



FIG. 4: Sample grayscale image (a) with regions of visibly 
distinct intensity and texture. Upper right (b) is the seg- 
mented image using the direct Potts method with Q = 16 
and St = 7.5. At lower left (c), the segmented image from 
(a) is filtered with a despeckle algorithm to reveal four dis- 
tinct texture categories. The categories are bright and smooth 
(white), bright and grainy (light gray), dark and grainy (dark 
gray), and dark and smooth (black). Lower right (d) is the 
despeckled segmentation result produced using the four-level 
histogram approach. The histogram method does not prop- 
erly distinguish regions of different textures, which demon- 
strates the problem in using the histogram approach for tex- 
ture analysis. 

1 or 3. Otherwise, (z, j) is assigned either texture 
category 2 or category 4). 

This algorithm guarentees that the total energy defined 
by Eq. [7] is minimized and processes in linear time since 
each pixel requires only a fixed number of computations. 
Only a single comparison and assignment required in 
Step 1, and a maximum of four comparisons (nearest 
neighbors) and assignment required in Step 2. Note that 
since the energy associated with each pixel is independent 
of the labeling of its neighbors, choosing an energy mini- 
mum for every pixel guarantees a total energy minimum 
(Eq. [7|) for the entire image. Finally, the segmented im- 
age is smoothed by applying the ImageMagick despeckle 
algorithm. This serves to eliminate much of the statisti- 



cal noise in the image. 



VI. RESULTS 

Three examples using the direct Potts segmentation 
will now be discussed using a photographic image, a med- 
ical image, and an acoustic seafloor image. Each segmen- 
tation result is compared to a histogram segmentation 
technique in which the image is partitioned into four uni- 
formly distributed intensity ranges. The four- level his- 
togram method is a very efficient, linear-time process but 
does not distinguish grainy regions from smooth regions. 
The photographic image in Fig. [5^ [29] contains depic- 
tions of water, sky, vegetation, and rocks. Using step 1 
of the direct Potts segmentation method as described in 
the previous section (with Q = 16 and St = 7.5), sky and 
rock are distinguished from water and vegetation accord- 
ing to pixel intensities. The sky and rock are represented 
as lighter shades, while water and vegetation are depicted 
as darker shades of gray in Figs. [SJd and [5t. To further 
separate the shades of gray, step 2 in the previous section 
effectively identifies regions of different texture (smooth 
or grainy). The full (steps one and two) direct Potts seg- 
mentation result is shown in Fig. [SJd. After filtering by 
using the despeckle algorithm, we see in Fig. [5t that the 
sky (white, bright /smooth texture) is distinguished from 
the rock regions (light gray, bright /grainy ) , and water 
(dark gray, dark/smooth) is distinguished from the vege- 
tation (black, dark/grainy). Segmenting the image by in- 
tensity using the four-level histogram produces the image 
in Fig. [5]i. The histogram approach does not distinguish 
sky from rock (both are light gray) and does compara- 
tively less well in distinguishing vegetation (lower left of 
image) from the background. A comparison of Figs. [St 
and [5]i underscores the importance of considering both 
intensity and texture in image segmentation. 

Figure [6t^[305 shows an MRI image of a brain with a 
tumor in the upper left quadrant. The segmented image 
using the direct Potts segmentation method with Q = 64 
and St = 31.5 is shown in Fig. [BJd. After filtering, we ob- 
tain the the image in Fig.[6t. Here, the white regions may 
be interpreted as possible tumorous tissue, light gray cor- 
responds to matter at the periphery of the tumor, dark 
gray represents cerebrum, and black is other brain mat- 
ter. Automated screening for tumors may be performed 
by finding white regions of significant size. Although, 
the four-level histogram method (Fig. [6li) also identifies 
a similar tumor region, the cerebrum is not distinguished 
from other brain matter. The affect of the chosen value 
for Q is demonstrated in Fig. [71 Larger values of Q re- 
sult in more of the image being categorized as textured 
(represented as light and dark gray). Finally, acoustic 
seafloor images may be analyzed in real-time with the 
direct Potts segmentation method as demonstrated in 
Fig. [SI Here sand ripples are identified (in green) as a 
grainy texture and is separated from the remainder of 
the image. Rapid seafloor characterization enables effec- 
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FIG. 5: Grayscale photo (a) showing water, sky, vegetation, 
and rocks. Segmentation of image using the direct Potts 
method (b) with Q = 16 and st = 7.5. Image is from the 
Berkeley Segmentation Data Set. (See 23.) The image at 
lower left (c) is a 4-level direct Potts segmentation of (a) that 
has been despeckled. Here the sky, rock, vegetation, and wa- 
ter are represented as white, light gray, dark gray, and black 
regions, respectively. The image at lower right (d) is the seg- 
mented and despeckled image using the four-level histogram 
method applied to (a). 




FIG. 6: MRI of a tumorous brain at upper left (a). A tumor 
is seen as the bright spot in the upper left region of the MRI. 
This image is from the Iowa Neuroradiology Library. Im- 
age at upper right (b) is the direct Potts segmentation with 
Q = 64, St = 31.5. Lower left (c) is the direct Potts segmen- 
tation that has been despeckled. The clearly visible white 
segment identifies the tumor. A tumor is clearly visible as 
the bright spot in the upper left region. The periphery of the 
tumor is shown in light gray, while cerebrum is represented 
as dark gray. Other brain matter is shown in black. For com- 
parison, lower right (d) is the despeckled four-level histogram 
segmentation, which identifies the tumor region, but cannot 
distinguish cerebrum from other brain matter. 



tive mine hunting/avoidance operations 0, 0, Hll, by 
indicating where seafloor mines are likely to be buried 
(mud), partially buried (sand), or unburied (rock). 



VII. CONCLUSIONS 

We have examined the direct application of the Potts 
model to grayscale image segmentation. The Potts model 
is applied to the grayscale image itself rather than its 
more common application to the segmented image. The 
result is a segmentation method which processes in lin- 



ear time instead of the exponential or power-law time 
required with conventional Potts segmentation methods. 
This increase in computational efficiency is vital for real- 
time, high data volume, and high resolution applica- 
tions. The method is analogous to minimizing the total 
energy in a two-dimensional magnetic spin glass where 
the energy is calculated according to four Potts mod- 
els. Each of the four Potts models corresponds to a dis- 
tinct type of magnetism, so the direct Potts segmenta- 
tion method effectively distinguishes regions in the spin 
glass that most closely adhere to the properties of four 
classes of magnetism, namely ferromagnetism, antiferro- 
magnetism, paramagnetism, and diamagnetism. In digi- 



(a) 



(b) 




FIG. 7: Despeckled segmentation images of Fig. [6^ using the 
direct Potts method are shown for different values for Q. Re- 
sults are shown for Q = 16 (a), Q = 32 (b), and Q = 128. 
In each case, the intensity threshold is st = (Q — l)/2. For 
higher values of Q, more of the image is categorized as tex- 
tured (light and dark gray). 




FIG. 8: Acoustic image of the seafloor formed using a sidescan 
sonar system, which creates an image going out on either side 
via an acoustic time series. The direct Potts segmentation 
with Q = 16 and st = 7.5 separates the region containing 
sand ripples (in green) from the rest of the image. The ripples 
are identified as a textured region and is labeled as sand since 
neither mud nor rock are likely to exhibit ripples. 



tal image segmentation, pixel intensities are substituted 
for the magnetic spins in the spin glass. The astute reader 
will note that no mention has been made of the physical 
basis for the coupling constants for the regions in the 
segmented image. The origins of the image pixels' cou- 
plings may, in fact, be complicated or unknown. Branch- 
ing of plants give a textured appearance in Fig. [5^, while 
the distribution of blood vessels and other brain material 
likely gives the textured appearance to Fig. [6^. Under- 
water sand ripples may arise from the fluid dynamics and 
saltation [26|, at the seafloor. The direct Potts seg- 
mentation has the advantages of being unsupervised (no 
training sets needed), very efficient (linear-time, only a 
few operation per pixel), and effective for texture iden- 
tification. For many image segmentation applications in 
which texture is important, this direct Potts method may 
be expected to produce better results than other linear- 
time methods. 



Acknowledgments 



This work was supported by the Naval Research Lab- 
oratory's Advanced Graduate Research Program. The 
author thanks the Tulane Center for Polymer Reac- 
tion Monitoring and Characterization (PolyRMC) for re- 
sources provided. A special thanks goes to Dr. Luca 
Celardo for helpful discussions. The views expressed in 
this article are those of the author and do not represent 
opinion or policy of the US Navy or Department of De- 
fense. 



[1] D. L. Pham, C. Xu, and J. L. Prince, "Current Methods 
in Medical Image Segmentation," Annu. Rev. Biomed. 
Eng. 2, 315-337 (2000). 

[2] Y. H. Yang, M. J. Buckley, S. Dudoit, and T. P. Speed, 
"Comparison of Methods for Image Analysis on cDNA 
Microarray Data," J. Comput. Graph. Stat. 11, 108-136 
(2002). 

[3] S. Peng, B. Urbane, L. Cruz, B. T. Hyman, and H. E. 



Stanley, "Neuron recognition by parallel Potts segmenta- 
tion," Proc. Natl. Acad. Sci. USA 100, 3847-3852 (2003). 

[4] V. Grau, A. U. J. Mewes, and M. Alcaniz "Improved 
watershed transform for medical image segmentation us- 
ing prior information," IEEE Trans. Med. Imaging 23, 
447-458 (2004). 

[5] S. Hadjidemetriou, C. Studholme, S. Mueller, M. Weiner, 
N. Schuff, "Restoration of MRI data for intensity non- 



7 



uniformities using local high order intensity statistics," 
Med. Image Anal. 13, 36-48 (2009). 
X. Descombes, M. Moctezuma, H. Maitre, and J. -P. 
Rudant, "Coastline detection by a Markovian segmenta- 
tion on SAR images," Signal Process. 55, 123-132 (1996). 
F. W. Bentrem, W. E. Avera, and J. Sample, "Estimating 
surface sediments using multibeam sonar," Sea Technol. 
47, 37-40 (2006). 

T. Asano, D. Z. Chen, N. Katoh, and T. Tokuyama, "Effi- 
cient Algorithms for Optimization-based Image Segmen- 
tation," Int. J. Comput. Geom. Ap. 11, 145-166 (2001). 

E. Ising, "Report on the theory of ferromagnetism," Z. 
Phys. 31, 253-258 (1925). 

R. B. Potts, "Some Generalized Order-Disorder Transfor- 
mations," Proc. Camb. Philos. Soc. 48, 106-109 (1952). 

F. W. Bentrem, "4-Level Spin-Glass Image Segmentation 
with Energy Minimization," Provisional Patent Applica- 
tion, Navy Case No. 99,755 (2009). 

K. Tanaka, "Statistical- mechanical approach to image 
processing," J. Phys. A: Math. Gen. 35, R81-R150 
(2002). 

J. P. Neirotti, S. M. Kurcbart, and N. Caticha, Phys. 
Rev. E 68, 031911 (2003). 

M. Blatt, S. Wiseman, and E. Domany, "Superparamag- 
netic Clustering of Data," Phys. Rev. Lett. 76, 3251- 
3254 (1996); "Data clustering using a model granular 
magnet," Neural Computation 9, 1805-1842 (1997). 
S. Wiseman, M. Blatt, and E. Domany, "Superparamag- 
netic clustering of data," Phys. Rev. E 57, 3767 (1998). 
K. Tanaka, H. Shouno, M. Okadak, and D. M. Tittering- 
ton, "Accuracy of the Bet he approximation for hyperpa- 
rameter estimation in probabilistic image processing," J. 
Phys. A: Math. Gen. 37, 8675-8695 (2004). 
E. Sharon, A. Brandt, and R. Basriy, "Fast Multiscale 
Image Segmentation," Comp. Vision Pattern Recogni- 
tion, Proc. IEEE Conf. 1, 70-77 (2000). 
P. F. Felzenszwalb and D. P. Huttenlocher, "Efficient 
Graph-Based Image Segmentation," Int. J. Comp. Vision 



59, 167-181 (2004). 
[19] A. X. Falcao, P. A. V. Miranda, and A. Rocha, "A Linear- 
Time Approach for Image Segmentation Using Graph- 
Cut Measures," Lec. Notes Comp. Sci. 4179/2006, 138- 
149 (2006). 

[20] D. Chandler, Introduction to Modern Statistical Mechan- 
ics (Oxford University Press, New York, 1987). 

[21] J. C. Lee Thermal Physics: Entropy and Free Energies 
(World Scientific Publishing Company, Singapore, 2002). 

[22] L. Onsager "Crystal statistics. I. A two-dimensional 
model with an order-disorder transition." Phys. Rev. 65, 
117-149 (1944). 

[23] D. Martin, C. Fowlkes, D. Tal, and J. Malik, A Database 
of Human Segmented Natural Images and its Application 
to Evaluating Segmentation Algorithms and Measuring 
Ecological Statistics, in Proceedings of the 8th Interna- 
tional Conference on Computer Vision, 2001 (2001), Vol. 
2, p. 416-423. 

[24] F. W. Bentrem, J. Sample, M. T. Kalcic, M. E. Duncan, 
High-frequency acoustic sediment classification in shal- 
low water, in Proceedings of Oceans 2002 (MTS/IEEE, 
Biloxi, 2002), Vol. 1, p. 7-11. 

[25] F. W. Bentrem, J. T. Sample, and M. M. Harris, Compu- 
tational Solutions for Today's Navy. Scientific Computing 
25, 30-32 (2008). 

[26] R. A. Bagnold The Physics of Blown Sand and Desert 
Dunes (Methuen, London, 1941). 

[27] L. Kang and L. Guo, Numerical simulation of aeolian 
sand ripples. Phys. Lett. A 330, 198-202 (2004). 

[28] ImageMagick Studio LLC, Landenberg, PA, USA 

[29] From the Berkeley Segmentation Data Set 

( http : / / www.eecs.berkeley.edu/Research/Projects /CS/ ~ vision / grouj 
Used with permission. 

[30] From the Iowa Neuroradiology Library 



(http:/ /www. uiowa.edu/c064s01/index.html ). 
with permission. 



Used 



